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ABSTRACT 

This  study  was  undertaken  to  investigate  the  effects 
of  varying  both  the  secondary  Mach  number  and  angle  of 
inclination  relative  to  the  primary  stream  on  the  flow 
field  generated  by  the  interaction  of  a  secondary  jet  with 
a  supersonic  mainstream.   The  experimental  portions  of 
this  investigation  were  conducted  at  a  primary  Mach  number 
of  2.80  in  the  Naval  Postgraduate  School  Supersonic  Wind 
Tunnel.   Data  are  presented  and  compared  with  various 
theories.   This  presentation  includes  correlation  of  the 
penetration  height  of  the  secondary  flow,  bow-shock  shape 
by  non-dimensionalization  with  respect  to  various  param- 
eters and  correlation  of  interaction  side  force  amplifi- 
cation factor. 
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I .   INTRODUCTION 

For  a  period  of  time  spanning  at  least  the  last  decade, 
considerable  interest  has  been  displayed  in  the  interaction 
of  secondary  gaseous  flows  injected  into  a  primary  super- 
sonic stream.   Several  theoretical  models  exist  which 
reflect  the  physical  trends  expected  for  varying  sets 
of  injection  parameters  (such  as  injection  pressure, 
temperature,  mass  flow  rate,  Mach  number  and  angle  of  incli- 
nation to  the  primary  flow,  et  cetera) .   The  body  of  experi- 
mental investigation  has,  for  the  most  part,  concerned 
itself  with  sonic  injection  normal  to  the  primary  stream. 

The  objective  of  this  work  is  to  study  the  phenomena 
associated  with  changes  of  both  injection  Mach  number  and 
inclination  angle  in  order  to  compare  results  with  the 
analytical  models  presently  available.   The  degree  to  which 
correlation  of  experimental  and  theoretical  results  are 
obtained  and  the  empirical  relationships  observed  will  be 
discussed. 

This  study  was  performed  at  the  Naval  Postgraduate 
School  Supersonic  Blow-Down  Wind  Tunnel  operating  at  a 
primary  Mach  number  of  2.80.   A  flat  plate  model  incor- 
porating an  interchangeable  nozzle  set  and  instrumented 
for  the  collection  of  plate  surface  pressure  distribution 
was  constructed  and  used  in  this  investigation.   Commer- 
cially available  dry  nitrogen,  helium  and  argon  were 
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injected   into  the  main  air  flow  at  various  stagnation 
pressures  and  mass  flow  rates  through  nozzles  of  varying 
Mach  number  and  angle  of  inclination  to  the  primary 
stream. 
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II.   EXPERIMENTAL  APPARATUS 

A.   THE  FLAT  PLATE  MODEL 

The  pylon  mounted  flat  plate  model  shown  in  Figure  1 
(made  of  steel  and  flash  plated  with  chromium)  was  mounted 
in  the  test  section  of  the  Naval  Postgraduate  School 
Supersonic  Wind  Tunnel.   The  tubing  which  leads  to  the 
forty-seven  pressure  taps  arranged  on  the  upper  surface 
of  the  flat  plate  passes  through  the  model  pylon  and  the 
tunnel  floor.   The  pressure  tap  locations  are  depicted 
in  Figure  2.   Access  for  mating  the  nozzles  with  the 
injectant  gas  supply  lines  is  also  made  available  through 
the  model  pylon  and  tunnel  floor.   A  representative  set  of 
drawings  of  three  interchangeable  injection  nozzles  is 
provided  in  Figure  3  showing,  from  left  to  right,  a  normal 
sonic  nozzle,  a  normal  supersonic  nozzle  and  an  inclined 
supersonic  nozzle.   All  nozzles  have  a  throat  diameter  of 
0.1000  ±  0.001  and  convergence  and  divergence  half-angles 
of  10°.   Nozzles  were  constructed  for  injection  at  Mach 
numbers  of  1.00,  1.73  and  2.13.   For  each  Mach  number 
nozzles  were  constructed  at  nominal  inclination  angles  e 
of  90°,  80°,  70°,  60°  and  50°.   Epsilon  is  taken  to  be  the 
angle  between  the  nozzle  axis  and  the  direction  of  the 
primary  stream  (or  the  surface  of  the  flat  plate)  in  a 
sense  such  that  for  injection  not  normal  to  the  primary 
flow  an  upstream  component  of  the  secondary  flow  exists. 
The  nozzles  are  inserted  from  above  the  model.   Alignment 
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and  a  positive  in-place  lock  for  securing  the  nozzles  in 
position  is  provided  by  a  small  key  bolt  passing  through 
the  pylon  from  below 

B.  THE  WIND  TUNNEL 

This  investigation  was  conducted  in  the  Naval  Post- 
graduate School  Blow-Down  Wind  Tunnel.   The  facility  has  a 
test  section  4  inches  by  4  inches  in  cross-section  by  6 
inches  in  length  and  provides  an  approximate  run  time  of 
5  minutes  at  a  Mach  number  of  2.80  and  a  stagnation  pressure 
of  50  psia.   The  lower  (tunnel  floor)  aluminum  block  was 
replaced  by  an  accurate  phenolic  replica  of  the  original 
lower  block  in  order  to  facilitate  mounting  the  model 
with  its  required  tunnel  floor  holes  while  preserving  the 
geometry  of  the  original  block  intact. 

C.  THE  EXPERIMENTAL  SET-UP 

Figure  4  is  a  schematic  representation  of  the  equipment 
arrangement.   The  calibrating  apparatus  is  not  shown. 

Injectant  gas,  manifolded  from  banks  of  four  commercial 
high  pressure  cylinders,  was  supplied  through  a  throttle 
valve  into  a  sharp -edged  orifice  flowmeter.   The  flowmeter 
upstream  pressure  was  measured  in  psia  by  a  Heise  gauge 
(Serial  No.  H  2946).   Differential  pressure  across  the 
orifice  was  sensed  by  a  Statham  differential  pressure 
transducer  (Serial  No.  11404)  calibrated  in  inches  of 
mercury  and  capable  of  measuring  a  ±  10  psi  differential 
at  a  working  pressure  of  500  psig  maximum.   The  stagnation 
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temperature  of  the  flow  was  measured  at  the  flowmeter  by 
means  of  a  chromel-alumel  thermocouple.   The  injectant  gas 
then  proceeded  to  an  expansion  chamber  mounted  adjacent  to 
the  lower  wall  of  the  wind  tunnel  in  which  the  flow  was 
effectively  stagnated.   At  this  station  the  total  pressure 
of  the  secondary  flow  was  measured  by  a  Heise  pressure 
gauge.   The  total  temperature  of  the  injectant  measured 
at  the  expansion  chamber  agreed  closely  with  that  measured 
downstream  of  the  flowmeter  (within  5°  F) .   The  stagnation 
temperature  at  the  flowmeter  was  taken  to  be  the  injection 
total  temperature.   The  total  temperature  and  total 
pressure  of  the  primary  flow  was  measured  in  the  wind  tunnel 
plenum  chamber. 

For  side  force  determination,  the  static  pressures  were 
measured  at  the  pressure  taps  on  the  flat  plate  model  by 
a  Scanivalve  (Model  48J4  -  1065  Scanivalve  with  a  0-10  psia 
transducer)  calibrated  in  inches  of  mercury  vacuum.   The 
measurements  were  recorded  on  a  Honeywell  2106  twelve 
channel  Visicorder  as  was  the  pressure  differential  across 
the  flowmeter  orifice  mentioned  above. 

Schlieren  photographs  on  Polaroid  4x5  type  55  P/N 
positive-negative  film  were  taken  of  each  run  using  an  off- 
set Schlieren  system  with  collimating  mirrors. 

Photographs  of  the  experimental  equipment  are  provided 
in  Figures  5  through  9. 
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III.   THEORETICAL  ANALYSIS 

From  the  many  phenomena  associated  with  secondary 
injection,  it  was  chosen  in  this  paper  to  deal  with  the 
characteristic  radius  which  may  be  computed  from  the  blast 
wave  analogy  model,  the  penetration  height  of  the  injected 
gas,  the  shape  of  the  interaction  or  bow  shock  and  the 
resultant  normal  side  force  developed  on  a  flat  plate. 

A.   PENETRATION  HEIGHT 

1 .   The  Blunt  Body  Model  of  Zukoski  and  Spaid 

The  similarity  of  the  shock  shapes  produced  by 
secondary  injection  into  supersonic  main  flows  and  the  shock 
patterns  produced  by  blunt  axisymmetric  bodies  in  supersonic 
flow  prompted  Zukoski  and  Spaid  [Ref.  1]  to  develop  a  blunt 
body  analogy  to  the  secondary  injection  phenomena.   Assuming 
the  proper  choice  of  a  blunt  body,  the  characteristic 
dimension  of  the  nose  (a  length  analogous  to  the  penetration 
height  of  the  secondary  injectant)  may  be  computed  by  equat- 
ing the  drag  of  the  required  blunt  body  to  the  change  in 
momentum  flux  of  the  secondary  injectant.   These  investi- 
gators assumed: 

a.  Sonic  injection  into  a  uniform  supersonic  flow 
with  no  boundary  layer 

b.  No  mixing  between  primary  and  secondary  flows 

c.  The  blunt  body  to  be  a  quarter  sphere  with  a 
semi-cylindrical  afterbody 
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d.   The  injectant  gases  remain  within  the  confines 
described  in  item  C 
By  assuming  Newtonian  flow  about  a  body  in  the  form  of  a 
quarter-sphere  with  a  half-cylinder  afterbody,  Zukoski  and 
Spaid  obtained  the  force  balance 


F  =  F  (1) 

x(Injection)     x(Blunt  Body) 


or 


m  V.2  =  (tt/8)  M2  Y  C*h2  P  (2) 

s  jx     x     °°  'p  p     °° 


V.   is  the  streamwise  velocity  of  the  injectant  and  is 
determined  by  assuming  isentropic  expansion  of  the  injectant 
from  the  secondary  stagnation  pressure  to  free-stream 

pressure.   C*   is  the  Newtonian  pressure  coefficient  behind 

*  P 

a  normal  shock  and  is  evaluated  at  the  nose  of  the  analogous 
blunt  body. 

Solving  Equation  (1)  for  the  radius  of  the  blunt 
body  (the  effective  penetration  height  h)  and  applying  a 
nozzle  coefficient  of  discharge  correction  to  permit  the 
use  of  physical  nozzle  dimensions  in  these  calculations, 
the  following  expression  was  obtained: 

h/d/5  =  w   (  >p    )  { (7-=T)  ty^) 

00   P   Pt  P      s     '  s 

PDt('Ys-1)/(Ys+1)  h  (3) 

•  [l-(p^)  s         s     ]}  (J) 

St 
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2 .   The  Blunt  Body  Model  of  Cassel,  Davis  and  Engh 

In  a  manner  somewhat  similar  to  the  method  of  Zukcski 
and  Spaid,  work  by  Cassel,  Davis  and  Engh  [Ref.  2]  employed 
an  equivalent  blunt  body  analogy  in  the  development  of  an 
expression  of  effective  penetration  height  of  a  secondary 
injectant.   The  obstacle  envisioned  in  their  model  is  a 
cylinder  which  is  tangent  to  a  flat  plate,  capped  on  the 
upstream  end  by  a  half-sphere  and  oriented  in  the  streamwise 
direction  of  the  primary  flow.   Assuming  that  the  momentum 
of  the  secondary  gases  is  parallel  to  the  plate  at  the 
sonic  line  in  the  plume  of  the  secondary  jet,  these 
investigators  also  equate  the  change  in  the  streamwise 
momentum  of  the  injectant  with  the  drag  of  the  equivalent 
body.   The  following  equation  results: 

C  A*qx  -  (P*-P.)  A*  =  m  [a*+V   sine(90°-e)]         (4) 

where  A*  is  the  projection  of  the  sonic  area  on  a  plane 
normal  to  the  wall,  P*  is  the  pressure  at  sonic  conditions, 
A  is  secondary  throat  area,  a*  the  speed  of  sound  at  sonic 
conditons  and  V   is  the  secondary  nozzle  exit  velocity. 


Continuity  implies 


2  .v<v1)p.A 


m  =  p*A*a*  =  p. A, a.  =  y  ( — -T) 

M         Kt  t  t    's  y  +1  a. 

s  t 

for  a  perfect  secondary  gas  at  constant  total  temperature. 
Therefore,  assuming  that  the  frontal  area  of  the  injectant 
stream  is  effectively  circular, 
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_L_  =  Si   =  -±   =  (_t)  (5) 

Pst    Pt    A*    lhj  ^ 


Substituting  into  Equation  (5)  the  two  previous 
equations  yields 

a*     h  2     ?   Y  /(Y  -DP  V 

k'^"^       "     (s#p7>{Ys'1+srin<90°-»)+1> 

In  Ref.  3,  the  vacuum  thrust  coefficient  of  a  nozzle  is  de- 
fined as 

Cf   "  'V1'^'  <6) 

o  '  s 

combining  the  above  equations,  Cassel,  Davis  and  Engh  obtain 

ql  Cd+(2/ysMoo2)   Ys  +  1     at  Ys  +  1 

where  the  last  two  terms  are  approximately  one  for  normal 
sonic  injection.   This  reduces  the  expression  for  sonic 
injection  to 


A*  =  F.  /q. 
jo  ^1 


or 


h  A*  ^       in   ^ 

dt     At       qxAt 


Crist,  Sherman  and  Glass  [Ref.  4]  obtained  a  good  fit  between 
analytical  and  experimental  results  with  the  expression 
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1 


hD/dt  =0.71  (Pst/PD) 


or,  from  Equation  (5) 


hD/h*  =0.71  (P*/PD) 


which  led  to  the  expression  for  e  =  90°. 


3 .   Penetration  Height  as  a  Function  of  Secondary 
Mach  Number 

Schetz  and  Billig  [Ref.  5]  propose  a  theoretical 
relationship  between  the  values  of  h  for  varying  M  which 
is  based  on  the  ratio  of  dynamic  pressure  at  initial  injec- 
tion conditions  to  that  at  mainstream  conditions.   For 
Y   =1.4  the  following  relationship  between  penetration 
height  at  M  =1.0  and  that  for  M   greater  than  1.0  was 
obtained 

.      1.2M  2 
h         s 


hl    1+0.  2  M  2 
s 


wh 


ere  h,  is  the  penetration  height  at  sonic  injection. 

B.   SHOCK  SHAPE  AND  SIDE  FORCE  PREDICTION 
1.   Blast  Wave  Analogy 

Hayes  [Ref.  6]  and  Lees  and  Kubota  [Ref.  7]  proposed 
a  similarity  between  hypersonic  flow  about  an  unyawed, 
axisymmetric  body  and  the  detonation  of  a  line  charge.   In 
general,  this  similarity  is  based  upon  correlating  the  time 
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dependent  pressure,  density  and  velocity  fields  which  are 
generated  by  the  explosion  of  a  line  charge  with  the 
similar  field,  generated  as  a  function  of  axial  distance 
over  the  velocity  of  the  undisturbed  flow,  about  a  body  in 
a  steady  hypersonic  flow.   In  their  blast  wave  analogy  j 

approaches  to  the  solution  of  the  secondary  injection  prob- 
lem, Broadwell  [Ref.  8]  and  Dahm  [Ref.  9]  have  employed 
this  similitude  concept  and  Sakuri ' s  [Refs.  10  and  11] 
solutions  to  the  flow  fields  behind  blast  waves  resulting 
from  the  explosion  of  line  charges  of  uniform  energy  per 
unit  length.   Two  basic  assumptions  are  common  to  the 
development  of  both  Broadwell  and  Dahm.   The  first  is  that 
the  injectant  gas  profile  is  equivalent  to  a  blunt  nosed 
body  and  the  second  that  the  pressure  field  between  the 
shock  and  a  blunt  nosed  body  is  correctly  predicted  by 
the  blast  wave  theory.   The  first  assumption  has  been 
justified  by  the  observed  similarity  of  the  pressure 
fields  about  blunt  bodies  to  those  generated  by  flat 
plate  injection.   The  second  is  less  tenable  in  that  it 
does  not  hold  strictly  true  where  hypersonic  flow  assump- 
tions do  not  apply,  i.e.,  near  the  nose  of  the  blunt  body 
or  equivalent  and  at  mainstream  Mach  numbers  not  in  the 
hypersonic  range.   Both  Dahm  and  Broadwell  assume  also  that 
at  some  point  not  far  down  stream  the  injectant  gas  has 
been  accelerated  to  the  mainstream  velocity. 

The  energy  per  unit  length  of  line  charge  in  the 
blast  wave  theoiy  may  be  considered  analogous  to  the  energy 
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added  to  the  flow  by  the  injectant.   Broadwell's  develop- 
ment considers  this  energy  term  to  be  composed  of  the 
secondary  momentum  change  resulting  from  acceleration  of 
the  injectant  to  the  mainstream  velocity  plus  the  energy 
generated  by  the  volume  addition  associated  with  injection. 
Dahm  derived  this  analogous  energy  term  directly  from 
thermodynamic  considerations  with  the  assumption  of  negli- 
gible mixing  of  primary  and  secondary  streams. 

In  the  blast  wave  theory,  a  full  360°  cylindrical 
shock  wave  is  considered  while  in  the  case  of  secondary 
injection,  where  the  flat  plate  is  a  plane  of  symmetry, 
only  a  half-cylinder  shock  wave  is  formed.   Therefore,  the 
analogous  energy  term  becomes 

E  =  2E' 

where  E  is  the  energy  per  unit  length  required  for  blast 
wave  formation. 

Sakuri  relates  a  characteristic  radius,  R  ,  to  the 
blast  wave  energy  and  primary  static  pressure  by: 

R  =  (E/2ttP  )%  (10) 

o         °° 

The  energy  term  as  derived  by  both  Broadwell  and  Dahm  is  of 
the  general  form 

E  =  2m  V  a).  (11) 

s  p  2 

Broadwell,  by  use  of  previously  mentioned  assumptions, 
ascribed  to  co^  the  expression, 
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V1   2 

"2  =  1+    (,  'J2){^7)  (12) 

(Yp-DM^     s  pt 

whereas,  from  thermodynamic  considerations,  Dahm  developed 

the    form 

Y  -1      2 

1  1  s       /  2       °°    .  ,    p    stN 

w2  =  I     ~~^ — r^T2  +  y-^r  ( —2 — >  (5mt:) 

Vv    «  Yp  °°         s 


Y  itl  T        ig  M  y-l9  y   -1      0   % 

+  'tWt1    5r  (1  +  -ErM-/1+-VMs  '    cose   (13> 

'p   s   pt  p 


where  e  is  the  angle  between  the  primary  flow  direction  and 
the  injection  axis. 

For  normal  injection,  with  which  Broadwell  was 
exclusively  concerned,  the  expressions  in  Equations  (12)  a 
and  (13)  have  been  noted  to  be  quite  similar  in  form 
although  their  origins  differ  considerably  in  concept. 
Henceforth,  oo^  will  be  used  here  exclusively  as  expressed 
by  Dahm  in  Equation  (13).   Therefore, 


u>~  =  E/2m  V 
2     '   s  p 


and   Equation    (10)    becomes 

k 

R     =    (<d_m  V  /ttP    )  (14) 

o  2    s   p        °° 

Taking 
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R  =  R/R 

'  o 

x  =  x/R 

'    o 

t  =  x/V 

P 

where  R  is  the  radius  of  the  interaction  shock  or  bow  shock 
at  x.   The  shock  radius  predicted  by  blast  wave  theory  is 
given  by  Sakuri  as 

a    ,  2  a   2        a   4 

(r^)  (4r)   =  J   [1  +  A   (-H  +  ^  L  (H   +  ...]    (15) 
Vj   R       o  Vj     2      Vj 

where  J   is  the  first  order  blast  wave  constant.  A.  are 
o  '   k 

higher  order  blast  wave  constants,  V  is  the  shock  velocity, 
and  a   is  sonic  velocity  in  the  undisturbed  mainstream. 
Solution  of  Equation  (15)  yields 

R  =  [2aR  /J  *]% 
o   o 

or,  for  the  first  order  theory, 

R  =  [2/Jo^]  2  (x/MJ3* 
and,  for  the  second  order  theory, 


2  ,  h    ,  2T  v  h    r ,     b   ,  x  ,  .,  ^ 


R=  ^'  (^)2  [1  +  7T^  'fr'1  (16> 

J        oo         2J  2    °° 
o  o 

The  theoretical  bow  shock  shape  for  any  given  set  of  injec- 
tion parameters  can  thus  be  determined. 

In  Broadwell's  development,  the  interaction  force 
is  described  by 
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F      =    2    fL    fR    (P-P    )    drdx 
w  o      o 


where 


P   =   P    [P    ,J    ,M    ,x,g(n) ] 

oo  '       o  °° 


which  yields 


F   =  F  [E,J  ,L,g(n)] 


w 


o 


where  g(n)  is  a  non-dimensional  pressure  ratio. 

In  a  similar  fashion,  Dahm  developed  the  expression 

for  second  order  theory  interaction  side  force 

3/2  1 

F.  12oo„         M   kn    ^  k      m   V      3/4 


(■ 


PA,  n       3/2  vd    '         P    A. 

°°        t  2    -r  t  °°        t 

TT     J 

o 
This  expression  applies  only  at  small  values  of  the  blast 
wave  analogy  parameter  s/R  M^  (less  than  0.0654  according 
to  Cassel,  Davis  and  Engh) ,   For  larger  values  of  this 
parameter,  i.e.,  for  lower  mass  flow  rates,  the  following 
relationship  was  developed  in  Ref.  9  for  circular  port 
injection 


F.  =  (2/7T)co„m  V  M  z  (17) 

l  2  s  p  °° 


where 


kl 
z  =  


t  (i  +  — °2  >  (4-)  /+£-  + 

2J  *    h%         knbJ/2  bh     V         b 
o  1 


,     J  k         ±  ±2 

1      o  ^ 


I 


n  (4-  +  ,  1+t-  )  ] 


Vs    kxb3/2'      b1"2 
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in  which  b  =  -J  A,,  and  <j>  =  bR  when  R  is  evaluated  at  s. 

Reference  9  evaluates  k,  =  0.4973  and  k~  =  -1.103  for 

Y   =  1.4. 
'P 

2 .   Upstream  Model 

Consideration  of  the  secondary  injectant  as  being 
analogous  to  a  blunt  nosed  obstacle  in  the  primary  stream 
lead  to  the  two-dimensional  upstream  model  of  Vinson, 
Amick  and  Liepman  [Ref.  12].   Here  the  weak  shock  formed 
by  boundary  layer  separation  in  the  high  pressure  region 
ahead  of  the  secondary  stream  or  the  analogous  blunt  body 
yields  the  equivalent  penetration  height  of  the  secondary 
injectant.   The  penetration  height  thus  determined,  the 
boundary  layer  detachment  angle,  the  mainstream  pressure 
and  the  pressure  behind  the  shock  in  turn  yield  the  magni- 
tude of  the  interaction  side  force  developed. 

Wu,  Chapkis  and  Mager  [Ref.  13]  extended  the  pre- 
vious model  to  three  dimensions  by  assuming  that  the 
separation  region  ahead  of  the  secondary  nozzle  is  conical. 
Their  model  of  the  interaction  phenomena  leads  to  the 
expression: 

2 
Y./y.    =   cos  a  tan  6  {cot  6  +  tan(oH-e')}  /2  y      /2/(Y  +D 

1    ~2  TT  S  S 

h 
,{(X  -  l)/[2Ys  +  (Ys-D  (j~  ~    1)]}         d8) 

where  a  is  the  primary  nozzle  half  angle  and  is  zero  for  the 
purposes  of  this  paper,  e'  =  90-e,  P  is  the  pressure  in  the 
area  of  the  detached  region  ahead  of  the  bow  shock,  6  the 
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boundary  layer  attachment  angle  and  6  the  associated  conical 
shock  angle. 

3 .   Downstream  Model 

Walker,  Stone  and  Shandor  [Ref.  14]  propose  a  model 
in  which  the  immediate  mixing  with  the  mainstream  gas  of 
a  low  flow  rate  injectant  in  a  constant  area  is  suggested. 
After  mixing,  the  gas  is  assumed  to  expand  isentropically 
to  the  mainstream  pressure  thus  inducing  a  side  force. 
The  following  expression  for  side  force  is  obtained: 

F.  cos    av      M  v    -1         n    -*g      y       h     m   T    ,  ^ 

_i= Ie_jl_  d+Iv  M   >    <-4r)  i-5^)  g     (i9) 


where 


Y~l         n         mmT  Y    ~1         oY         Y    _1 

2  -  ms      msTpt  2  -         Yp     V1 


T    "I 
m  m  T    ,  h       2y       h       1+  -% —  M 


I    A     P M    2    ^ 


in  l  4n  T    . J       S    +1  T  ' 

s  s   pt  's  'p 


U    +    (Y    _!)    M   2]    sin    (a+e) 


p  oo  M 
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IV.   DESCRIPTION  OF  THE  EXPERIMENT 

A.   PARAMETERS 

The  secondary  nozzles  constructed  for  use  in  this 
experiment  were  of  Mach  numbers  1.00,  1.73  and  2.13.   A 
description  of  the  calibration  of  the  supersonic  nozzles 
is  contained  in  Appendix  A.   The  throat  diameter  of  all 
nozzles  was  a  nominal  0.100  ±  0.001  inches.   Nominal 
inclination  angles  of  90°,  80°,  70° ,  60°  and  50°  were 
tested  at  each  injection  Mach  number. 

During  all  runs,  the  primary  flow  Mach  number  remained 
constant  at  M  =2.80.   The  primary  total  pressure  was  a 
constant  50.0  psia.   In  successive  closely  spaced  runs, 
the  primary  total  temperature  decreased  slightly  as  the 
tunnel  walls  were  cooled  by  the  flow.   This  temperature 
decrease  was  slight,  having  an  overall  variance  of  ±  15° F 
from  the  mean  total  primary  temperature  for  all  runs.   For 
purposes  of  the  interpretation  of  the  results  of  the 
theoretical  computations  involving  the  primary  flow 
parameters,  this  temperature  variation  was  considered 
negligible . 

The  secondary  total  pressure  was  increased  by  increments 
of  25.0  psia  from  25.0  to  250.0  psia  providing  a  series  of 
ten  secondary-to-primary  total  pressure  ratios  for  each  of 
the  15  nozzles  tested.   The  secondary  gas  reservoir  total 
temperatures  varied  only  slightly  due  to  minor  fluctuations 
in  the  ambient  temperature  over  the  period  of  the  testing. 
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Therefore,  the  total  temperature  of  the  secondary  flow  after 
it  was  throttled  to  the  desired  total  pressure  also  varied 
slightly.   This  effect  was  also  neglected  for  purposes  of 
interpretation . 

B.   EXPERIMENTAL  PROCEDURE 

The  test  runs  were  conducted  in  ascending  order  of 

M   and  at  a  given  secondary  Mach  number  in  descending  order 

of  e .   Ten  runs  were  conducted  at  each  nozzle  configuration 

at  the  previously  described  total  pressure  ratios  using 

nitrogen  as  the  secondary  gas.   Two  nozzle  configurations 

(M   =  1.0,  e  =  90°  and  M   =  1.73,  e  =  68°)  were  tested 
s  s 

similarly  with  argon  and  helium  as  the  secondary  injectant. 

Individual  runs  averaged  less  than  one  minute  in 
duration.   During  each  run  the  following  observations  were 
made : 

1.  The  static  pressure  at  each  of  the  47  flat  plate 
pressure  taps  was   measured  by  a  Scanivalve  and 
recorded  on  a  Visicorder. 

2.  The  differential  pressure  through  the  flowmeter 
was  measured  by  a  differential  pressure  trans- 
ducer and  recorded  on  a  Visicorder. 

3.  The  upstream  flowmeter  pressure  was  measured  by 
pressure  gauge. 

4.  The  injection  total  pressure  was  measured  by 
pressure  gauge , 
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5.  The  primary  total  pressure  was  measured  by  pressure 
gauge . 

6.  The  total  temperature  of  the  secondary  flow 
through  the  flowmeter  was  sensed  by  a  chromel- 
alumel  thermocouple  and  displayed  on  a  millivolt- 
meter. 

7.  Primary  total  temperature  was  sensed  by  a  chromel- 
alumel  thermocouple  and  displayed  on  a  millivolt- 
meter. 

8.  A  Schlieren  photograph  was  taken  of  the  test 
section. 

As  described  earlier,  the  total  temperature  of  the  secondary 
flow  at  the  flowmeter  was  taken  to  be  the  injection  total 
temperature. 

Tare  pressures  were  observed  at  the  plate  surface  with 
the  tunnel  operating  but  without  secondary  injection. 
These  were  then  subtracted  from  the  respective  pressures 
measured  during  injection  to  determine  the  interaction 
pressure  rise.   This  pressure  rise  was  then  integrated 
over  the  plate  surface  for  determination  of  the  interaction 
side  force.   Tare  observations  were  made  prior  to  each  of 
the  series  of  10  runs  for  a  given  nozzle  configuration. 
In  all  cases,  tare  pressures  were  observed  with  a  nozzle 
for  normal  injection  installed. 

Penetration  heights  and  bow  shock  coordinates  were 
measured  by  projecting  the  negatives  of  Schlieren  photo- 
graphs onto  a  measuring  grid.   The  grid  was  constructed  to 
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twice  the  actual  scale  of  the  model  and  reference  points 
were  provided  to  insure  that  the  grid  was  properly 
oriented  with  respect  to  the  projection  of  the  negative. 

In  order  to  determine  the  effect  of  the  protrusion  of 
the  downstream  side  of  the  nozzles  at  E  <  90° ,  a  duplicate 
sonic  normal  nozzle  was  constructed.   A  dummy  obstruction 
representing  the  protrusion  of  an  inclined  nozzle  was  then 
inserted  just  downstream  of  the  normal  sonic  nozzle  orifice 
This  nozzle  was  tested  in  order  to  determine  the  effect 
of  the  obstruction  on  the  pressure  distribution  over  the 
plate  and  on  penetration  height. 

A  normal  sonic  nozzle  was  constructed  with  a  cylin- 
drical constant  area  throat  section  0.01  inches  in  length 
and  tested  for  comparison  with  the  sharp-edged-throat 
nozzles  used  in  this  experiment 
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V.   EXPERIMENTAL  RESULTS 

A.   PENETRATION  HEIGHT 

Zukoski  and  Spaid  have  stated  that  the  height  of 
penetration  of  the  secondary  gases  into  the  main  flow  is 
a  dimension  characteristic  of  the  scale  of  the  interaction 
disturbance.   This  parameter,  non-dimensionalized  by  the 
nozzle  throat  diameter  and  the  square  root  of  the  nozzle 
discharge  coefficient  (i.e.,  by  the  effective  nozzle 
throat  diameter)  has  been  utilized  by  previous  investiga- 
tors [Refs.  1,  15  and  16]  to  correlate  injection  penetration 
height  with  secondary-to-primary  total  pressure  ratios  for 
normal,  sonic  injection. 

In  this  experiment,  the  discharge  coefficient  was  found 
by  dividing  the  actual  flow  rates  through  the  nozzles  by 
the  corresponding  mass  flow  rate  calculated  for  ideal  gas 
assuming  isentropic  flow.   The  actual  flow  was  measured  by 
means  of  a  sharp-edged  orifice  flowmeter  placed  in  the  injec- 
tant  supply  line.   The  flowmeter  was  constructed  and  the 
flow  through  it  evaluated  in  accordance  with  Ref.  17. 
Discharge  coefficients  thus  obtained  are  shown  versus  the 
inverse  square  root  of  nozzle  throat  Reynolds  number  in 
Figures  10,  11,  and  12.   These  results  are  in  general 
agreement  with  the  discharge  coefficients  of  similar  nozzles 
determined  by  Amick  and  Hays  [Ref.  16]  and  by  Chrans  and 
Collins  [Ref.  15]. 
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As  described  in  detail  earlier,  measurements  of  the 
penetration  heights  were  taken  from  the  Schlieren  photo- 
graphs which  were  made  of  each  run.   This  dimension  was 
measured  along  a  normal  from  the  surface  of  the  flat  plate 
to  the  point  of  highest  penetration  of  the  secondary  flow. 

For  injection  at  a  given  Mach  number,  penetration  height 

correlations  in  non-dimensional  plots  of  (h/d/c)  sin  e  versus 

(P  ,/P  , )  sin  e  for  all  £  were  obtained.   The  results  are 
st   pt 

shown  in  Figures  13,  14  and  15.   For  the  special  case  of 

normal  (sin  e  =  1)  ,  sonic  injection  Figure  13  is  equivalent 

to  a  plot  of  h/d/c   as  a  function  of  P  ,/P   which,  as 
c  st   pt 

stated  earlier,  has  seen  previous  use.   The  present 
application  extends  the  generality  of  the  earlier  correlation. 
From  Figures  13,  14  and  15  it  can  be  seen  that 

p 
,  h  s  r  ,    st,.       -,  3 

( )  sin  e  =  [  (- — )  sin  £] 

d/c  pt 

where  the  exponent  6  varies  with  Mach  number.   In  Figure 
13,  3  was  found  to  be  approximately  equal  to  0.53.   This 
was  in  close  agreement  with  the  value  of  3  =  0.5  in  the 
theoretical  development  of  Zukoski  and  Spaid  for  normal 
sonic  injection.   As  can  also  be  seen  in  Figure  13,  the 
observed  data  agrees  closely  with  the  theoretical  pene- 
tration height  for  normal,  sonic  injection  of  Cassel,  Davis 
and  Engh.   The  Riemann  disc  height  as  a  function  of  pressure 
ratio  for  an  underexpanded  sonic  nozzle  exhausting  into 
still  air  as  determined  in  Ref.  18  is  also  shown  in  Figure 
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13.   This  curve  agrees  closely  with  both  the  observed  and 
the  theoretical  correlations. 

As  expressed  by  Schetz  and  Billig  in  Equation  (9),  the 
penetration  height  is  expected  to  increase  as  a  function 
of  Mach  number.   In  Figures  13,  14  and  15  this  is  borne  out 
in  a  qualitative  sense.   Figure  16  is  a  reproduction  of  the 
mean  lines  defined  by  the  data  in  the  preceding  three 
figures.   As  shown  in  this  figure,  the  penetration  height 
does  increase  with  increasing  secondary  Mach   number. 
However,  values  of  the  exponent  3  are  also  seen  to  increase. 
As  indicated  by  the  theoretical  plots  also  shown  in  Figure 
16,  this  was  not  predicted  by  Schetz  and  Billig.   At 
M  =1.73  the  observed  value  of  $  is  approximately  0.62  and 
at  M  =2.13  approximately  0.63. 

B.   SHOCK  SHAPE 

Dahm  and  Broadwell  suggest  that  the  characteristic  radius 

R   is  a  dimension  representative  of  the  scale  of  the  inter- 

o  r 

action.   Dahm  has  shown  theoretically  that  a  non-dimensional 
plot  of  x/R  =  x  versus  R/R  =  R  provides  correlation  of 
bow-shock  coordinates.   Experimental  verification  of  this 
correlation  for  normal  sonic  injection  was  obtained  by 
Chrans  and  Collins  with  variations  in  mass  flow  rate  and 
secondary-to-primary  total  temperature  ratio.   The  present 
study  has  shown  that  these  coordinates  are  not  sufficiently 
general.   Simultaneous  correlation  is  not  achieved  when  data 
is  collected  at  various  values  of  inclination  angle  e.   As 
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shown  for  nitrogen  injection  in  Figures  17  through  30/  it 
was,  however,  found  that  bow  shock  shape  correlated  in  a 
plot  of  x  versus  R  for  any  specific  E.   Correlation  with 
respect  to  all  values  of  epsilon  was  obtained  by  plotting 
shock  shapes  in  the  non-dimensional  coordinates  x  sin  e 
versus  R  sin  e.   For  normal  injection,  where  sin  e  =  1,  the 
coordinates  in  effect  revert  to  x  and  R.   Bow  shock  shapes 
are  correlated  in  this  manner  in  Figures  31,  32,  and  33  for 
nitrogen,  argon  and  helium  injected  at  a  secondary-to- 
primary  total  pressure  ratio  of  3.0.   It  can  be  seen  from 
the  previous  data  correlations  for  various  pressure  ratios 
in  Figures  17  through  30  that  the  results  here  are  repre- 
sentative of  all  pressure  ratios. 

The  data  scatter  in  Figures  31,  32,  and  33,  although 
coherent,  is  separated  to  a  slight  degree  with  respect  to 
secondary  Mach  number.   In  general,  the  data  representing 
injection  at  M  =2.13  lie  along  the  upper  bound  of  the 

scatter  band,  that  for  M  =  1.73  in  the  center  and  the  data 

s 

for  sonic  injection  along  the  lower  bound.   The  parameter 

R  ,  by  which  bow-shock  coordinates  were  non-dimensionalized 

o    ■* 

in  these  figures,  was  calculated  from  Equation  (14)  in 

which  R   is  a  function  of  u)„  .   As  can  be  seen  from  Equation 
o  2  ^ 

(13),  a)^  is  a  function  of  M   only  for  e  <  90°.   This  depen- 
dence on  M   is  sufficient  to  cause  shock  shapes  to  correlate 
very  closely  for  all  e  at  a  given  M  .   Correlation  at 
differing  values  of  M   is,  however,  not  so  precisely 
obtained.   In  the  previous  section,  h  was  suggested  as  a 
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characteristic  dimension  representative  of  the  scale 
of  the  disturbance.   In  Refs.  1  and  15  bow  shock  shape 
correlation  for  the  case  of  normal,  sonic  injection  was 
obtained  by  plotting  x/h  =  x'  versus  R/h  =  R1 .   Corre- 
lation of  shock  shapes  for  all  inclination  angles  was 
obtained  in  Figure  34  by  plotting  R1  sin  e  versus  x'  sin  e. 
For  normal  injection  (sin  e  =  1)  these  coordinates  are 
identical  with  the  coordinates  previously  used  for  correla- 
tion of  shock  shapes  formed  by  normal  sonic  injection.   It 
was  found  that  when  correlated  in  this  manner,  the  data 
for  various  secondary  Mach  numbers  was  randomly  distributed 
throughout  the  scatter  band  rather  than  forming  obvious  sub- 
groups by  M  within  it  as  was  the  case  when  shock  coor- 
dinates were  non-dimensionalized  by  R  . 

■*   o 

The  mean  lines  representative  of  the  shock  shape  cor- 
relation data  in  Figures  31,  32,  and  33  (for  nitrogen, 
argon  and  helium)  are  reproduced  for  comparison  in  Figure 
35.   The  dependence  of  shock  shape  on  molecular  weight 
previously  noted  by  Chrans  and  Collins  is  seen.   At  a 
given  value  of  x  sin  e  the  non-dimensional  radius,  R  sin  e 
increased  with  increasing  secondary  molecular  weight.   Blast 
wave  theory  does  not  indicate  such  a  dependence.   It  has 
been  suggested  [Ref.  15]  that  the  mixing  of  the  secondary 
gases  with  the  primary  stream  (assumed  to  be  negligible  in 
the  development  of  the  blast  wave  analogy)  may  account  for 
this  dependence. 
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When  plotted  in  terms  of  x1  sin  e  versus  R1  sin  e  ,  the 
same  molecular  weight  dependence  is  shown. 

C.   SIDE  FORCE 

An  experimental  examination  of  the  interaction  side 
force  developed  by  secondary  injection  was  also  conducted. 
Side  forces  on  the  flat  plate  were  calculated  from  the 
interaction  pressure  rise  on  the  plate  surface  as  detailed 
earlier.   The  interaction  side  force  was  calculated  for  all 
runs . 

In  this  study,  side  forces  are  presented  in  the  form 
of  a  non-dimensional  force  amplification  factor  A.   The 
amplification  factor  is  defined  as 
F.  +  F .  sin  e 


A 


1 1 


F  .  sin  e 
3 


which  relates  the  interaction  side  force  to  the  component 
of  secondary  nozzle  thrust  normal  to  the  plate  surface. 
Reproducibility  of  the  observed  values  of  A  were  consis- 
tently good.   Figure  36  (a)  shows  representative  data 
scatter  in  a  plot  of  A  versus  primary-to-secondary  total 
pressure  ratio. 

In  order  to  prevent  disruption  of  the  symmetry  of  the 
secondary  flow  by  the  physical  hardware  of  the  model,  all 
nozzles  were  constructed  with  the  upstream  edge  of  the 
nozzle  exit  planes  coincident  with  the  plate  surface.   As 
a  result,  the  downstream  edge  of  the  nozzles  protruded 
slightly  into  the  primary  flow.   It  was  necessary  to 
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determine  the  effect  of  this  protrusion.   A  sonic,  normal 
nozzle  was  constructed  with  an  obstruction  representing  an 
inclined  nozzle  installed  immediately  downstream  of  the 
nozzle  orifice.   This  nozzle  was  then  tested  for  comparison 
with  the  normal  sonic  nozzle  used  in  the  experiment.   The 
values  of  A  at  various  injection  pressure  ratios  thus 
obtained  agreed  within  the  normal  experimental  scatter 
encountered.   A  normal  sonic  nozzle  having  a  short  (0.01 
inches)  cylindrical,  constant  area  throat  (square  edged) 
was  also  tested  for  comparison  with  the  sharp  edged  nozzles 
used  in  this  experiment.   The  amplification  factors  were 
found  to  agree  within  the  bounds  of  anticipated  scatter. 
The  results  of  these  comparisons  are  shown  in  Figure  36  (b) 
Schlieren  photographs  showed  no  difference  in  either  bow 
shock  shape  or  penetration  height  in  these  comparisons  of 
nozzle  configuration. 

In  the  development  of  their  blunt  body  models,  Zukoski 
and  Spaid,  and  others  [Refs.  2  and  8]  suggest  that  the 
characteristic  length  is  a  function  of  the  change  in 
momentum  experienced  by  the  secondary  gases. 

As  has  been  shown  earlier,  penetration  height  increases 
as  epsilon  decreases  (i.e.,  as  the  upstream  component  of 
momentum  increases)  when  all  other  parameters  are  held  con- 
stant.  It  is,  therefore,  to  be  expected  that  as  the  scale 
of  the  disturbance  increases  the  interaction  side  force 
will  increase  accordingly.   In  Figure  37,  at  a  repre- 
sentative secondary-to-primary  total  pressure  ratio  of 
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3.0,  a  plot  of  interaction  side  force  versus  epsilon  is 

shown.   It  may  be  seen  that  a  decrease  in  epsilon  is 

accompanied  by  an  increase  in  observed  side  force.   This 

figure  also  shows  a  decreasing  trend  in  interaction  side 

force  as  M   increases.   The  interaction  side  force  for 
s 

sonic  injection  and  that  for  injection  at  M   =  1.73  at 
J  -J  s 

various  values  of  e    are  nearly  equal.   For  injection  at 

M  =2.13,  however,  a  decrease  in  observed  side  force 
s        '         ' 

occurs.   This  decrease  is  not  predicted  by  Dahm.   Reference 
5  suggests  that  the  gains  which  theoretically  appear  to  be 
associated  with  supersonic  injection  may  be  countered  by 
the  losses  involved  in  turning  a  supersonic  secondary  jet 
downstream.   The  turning  of  a  sonic  secondary  jet  is 
largely  shock  free  and  the  losses  are  thus  small  compared 
to  the  turning  of  a  supersonic  jet. 

In  Figures  38,  39,  and  40,  interaction  side  force  is 
plotted  in  terms  of  the  nondimensional  force  amplification 
factor  A  versus  secondary-to-primary  total  pressure  ratio. 
The  amplification  factor  A,  at  a  given  value  of  pressure 
ratio,  is  seen  to  increase  with  decreasing  epsilon.   Part 
of  this  increase  is  due  to  the  decreasing  value  of  the 
normal  component  of  jet  thrust  (F .  sin  e)     as   e  decreases 
and  part  is  due  to  the  increasing  value  of  F.  as  e  decreases 
which  is  shown  in  Figure  37. 

It  is  noted  in  Ref.  19  that  an  empirical  correction 
factor  of  0.51  when  applied  to  Dahm's  second  order  theory 
values  of  side  force  provides  a  semi-empirical  value  of 
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side  force  very  nearly  comparable  with  previous  experimental 
data.   The  shaded  areas  on  Figures  38,  39  and  40  represent 
the  region  of  the  value  of  A  computed  from  this  semi- 
empirical  correction  to  Dahm's  theory.   It  can  be  seen 
that  considerable  data  scatter  exists  in  the  area  of  low 
injection  pressure  ratios  and  that  observed  and  semi- 
empirical  values  differ  considerably.   Conversely,  at  the 
higher  pressure  ratios,  scatter  decreases  and  agreement 
with  the  semi-empirical  values  of  A  improves  considerably. 
Also,  as  M   increases,  a  comparison  of  Figures  38,  39  and 
40  shows  that  overall  data  scatter  decreases  and  overall 
agreement  with  the  semi-empirical  values  of  A  increases. 

In  Figures  41,  42  and  43,  the  data  shown  in  Figures 
38,  39  and  40  is  presented  in  a  more  general  manner.   The 
previously  observed  successes  in  correlation  of  parameters 
suggest,  as  e  is  varied,  that  all  features  of  the  flow 
scale  by  sin  e  when  properly  non-dimensionalized.   A  degree 
of  correlation  was  achieved  in  Figures  38,39,  and  40  in 
a  plot  of  A  sin  e   versus  m  V.  sin  e.   Again,  at  lower  values 

of  m  V.  sin  e(i.e.,  lower  injection  pressure  ratios)  data 

s  j  J  c 

scatter  is  pronounced  and  agreement  with  empirically 
corrected  theory  is  not  good.   As  m  V .  sin  e  increases, 
however,  data  scatter  decreases  and  agreement  with  the 
corrected  theory  continually  improves.   As  M  increases, 
comparison  of  Figures  41,  42  and  43  also  shows  considerable 
improvement  in  the  degree  of  correspondance  of  data  with 
the  curve  representing  the  semi-empirical  values  of  A. 


42 


The  degree  to  which  observed  values  of  side  force  agree 

with  the  side  force  actually  developed  is  a  function  of  the 

fineness  of  the  grid  spacing  of  the  pressure  taps.   It 

was  found  in  this  experiment  that  both  data  scatter  and 

disagreement  with  corrected  theory  were  pronounced  at  the 

lower  values  of  injection  pressure  ratio  (or  the  related 

lower  values  of  m  V.  sin  e).   In  this  area  of  operation, 

s  j  r 

the  characteristic  length  h  was  of  the  order  of  or  less 
than  the  average  distance  between  pressure  taps.   As  h 
increased  to  a  value  of  the  order  of  twice  the  pressure 
tap  spacing  or  greater  (i.e.,  for  values  of  (P   /P  ,)  sine 
greater  than  2.0)  data  scatter  decreased  and  agreement  with 
semi --empirically  predicted  values  of  A  increased.   Above 
values  of  (P  ,/P   )  sin  e    >  2.0  or  the  corresponding  values 
of  m  V.  sin  e  >  25,  the  data  in  Figures  38  through  43  is 
considered  reliable.   In  this  area,  the  data  in  Figures 
41,  42,  and  43  shows  a  consistant  degree  of  correlation 
which  agrees  well  with  expected  values. 

As  determined  by  Zukoski  and  Spaid,  the  surface  area 
ahead  of  the  injector  is  subjected  to  the  highest  inter- 
action pressures  developed  during  secondary  injection. 
It  was  observed  in  this  experiment  that  the  pressures 
developed  downstream  of  the  injector  were  generally  of 
little  significance  in  calculating  the  value  of  inter- 
action side  force.   The  model  of  the  interaction  proposed 
by  Wu,  Chapkis  and  Mager  depends  exclusively  upon  flow 
phenomena  upstream  of  the  nozzle  for  determination  of  the 
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interaction  side  force.   A  limited  investigation  of  this 
theory  was  undertaken.   The  results  of  the  comparison  of 
observed  data  and  that  calculated  according  to  Equation 
18  is  shown  in  Figure  44  in  a  plot  a  A  versus  P   /P   . 
The  calculated  values  were  determined  by  measuring  the 
boundary  layer  detachment  angle  6  from  the  Schlieren  photo- 
graphs of  the  respective  runs.   These  measurements  were 
accomplished  on  an  optical  comparitor.   The  shock  angle  6 
was  not  shown  on  the  photographs  with  sufficient  clarity 
to  permit  measurement.   Therefore,  the  value  of  9  corres- 
ponding to  a  cone  of  half-angle  6  was  taken  from  Ref.  20. 
The  pressure  behind  the  conical  shock  was  also  determined 
from  Ref.  20.   The  results  of  the  computations  thus 
achieved  are  shown  by  the  curve  in  Figure  44.   The  broken 
portion  of  the  curve  represents  areas  in  which  photographic 
quality  did  not  permit  accurate  measurement  of  6.   Observed 
and  calculated  values  agreed  closely. 

The  side  force  determined  by  the  downstream  theory  of 
Walker,  Stone  and  Shandor  was  evaluated  in  accordance  with 
Equation  19  for  normal  sonic  injection.   The  value  of  A 
at  low  values  of  pressure  ratio  was  found  to  be  of  the 
order  of  6,  in  considerable  disagreement  with  observed  and 
previously  computed  values. 

In  Figures  45  and  46,  values  of  A  are  plotted  against 
pressure  ratio  and  m  V.  sin  e,  respectively,  for  nitrogen, 
argon  and  helium  for  two  different  nozzle  configurations. 
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Within  the  previously  encountered  range  of  data  scatter, 
no  dependence  of  A  on  the  molecular  weight  of  the  injectant 
was  observed. 
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VI.   CONCLUSIONS 

Penetration  height  increases  both  with  an  increase  in 

M  and  with  a  decrease  in  e,    all  other  parameters  being 

held  constant.   Correlation  of  penetration  height  for  a 

given  value  of  M   is  obtained  in  a  plot  of  (h/d/c)  sin  e 

versus  (P   /P   )  sin  e.   For  normal  sonic  injection  these 

coordinates  become  the  familiar  h/d/c  versus  P  ./P  .  .   In 

st7  pt 

the  case  of  sonic  injection  ,  the  slope  of  the  correlation 
curve  is  approximately  0.5,  i.e., 

.  P  .         0.5 

h     •    ^  r  /  StN  , 

sine^t  (- — )  sin  e  J 

d/c         pt 

or  for  the  case  of  normal  sonic  injection 

0.5        0.5 
h/d/c"  *[Pst/Ppt]     a  (ms) 

This  is  in  close  agreement  with  the  theories  discussed. 
As  M  increases,  the  value  of  the  exponent  increases.   At 
M  =1.0  the  value  of  the  exponent  is  approximately  0.53, 
at  M  =1.73  approximately  0.62  and  at  M  =2.13  approx- 

S  o 

imately  0.63. 

Non-dimensional  bow  shock  shape  plotted  in  the  coordinate 
set  (x,  R)  remains  essentially  constant  with  varying  M  . 
The  non-dimensional  radius  at  a  given  value  of  x  (where  x 
is  measured  downstream  from  the  nozzle  centerline)  increases 
with  decreasing  e.   Bow  shock  shapes  correlate  in  non- 
dimensional  plots  of  x  sin  e  versus  R  sin  e  or  x1  sin  e 
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versus  R'  sin  e.   The  correlation  curves  are  coincident 
with  those  for  normal  sonic  injection  plotted  in  the  coor- 
dinates x  versus  R  or  x1  versus  R'. 

The  empirical  correction  factor  of  0.51  applied  by 
Dahm  to  the  values  of  interaction  side  force  calculated 
from  his  second  order  theory  was  qualitatively  substan- 
tiated.  Correlation  of  side  force  data  in  a  plot  of  A 

sin  e  versus  m  V.  sin  e  is  indicated  for  all e  and  M  . 
s  j  s 

In  order  to  obtain  valid  side  force  data  by  means  of  inte- 
grating interaction  pressure  rise  over  a  flat  plate,  the 
spacing  of  the  pressure  taps  should  be  on  the  order  of  or 
less  than  one-half  the  smallest  penetration  height  expected. 

Dahm's  theory  predicts  an  increase  in  the  magnitude  of 
side  force  with  decreasing  e,  all  other  parameters  remain- 
ing constant.   This  was  qualitatively  confirmed.   Side 
force  at  a  given  value  of  e   was  observed  to  decrease  with 
increasing  M  .   Dahm's  theory  predicts  no  change  with 
increasing  M   for  normal  injection  and  a  slight  increase 

in  side  force  with  increasing  M   as  e   decreases. 

r   s 

For  the  same  values  of  injectant  momentum,  the 
experimental  results  indicate,  within  the  degree  of  data 
scatter  expected,  that  side  force  does  not  vary  with  changes 
in  the  molecular  weight  of  the  secondary  gas.   At  the  nozzle 
throat,  i.e.,  at  choked  flow,  the  secondary  mass  flow  rate 
is  reduced  as  the  molecular  weight  of  the  injectant  de- 
creases.  It  is  beneficial,  therefore,  to  inject  gases  of 
low  molecular  weight  in  order  to  minimize  the  total  weight 
of  secondary  fluid  which  must  be  carried  by  the  vehicle. 
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APPENDIX  A 
SUPERSONIC  NOZZLE  CALIBRATION 

The  supersonic  nozzles  utilized  in  this  experiment  were 
constructed  with  10°  divergence  half-angles  on  either  side  of 
a  sharp  edged  throat.   The  diameters  of  the  throats  of  all 
nozzles  were  0.10  ±  .001  inches.   Throat  diameters  were  held 
to  this  tolerance  in  the  following  manner.   Two  carefully 
manufactured  rods  were  constructed,  one  having  a  diameter 
of  0.10  -  .001  =  0.099  inches  and  the  other  having  a  diameter 
of  0.10  +  .001  =  0.101  inches.   The  10°  half-angle  approach 
portion  of  the  nozzle  was  drilled  out.   The  portion  of  the 
nozzle  downstream  of  the  throat  was  then  drilled  (also  with 
the  10°  half-angle  bore)  until  the  smaller  diameter  rod  but 
not  the  larger  diameter  rod  was  free  to  pass  through  the 
nozzle  throat.   Thus,  the  throat  diameter  was  within  the 
diameter  difference  of  the  testing  rods,  or  0.10  ±  .001 
inches . 

The  nozzle  having  the  lowest  supersonic  Mach  number  was 

constructed  to  a  nominal  area  ratio  of  1.6900  corresponding 

to  M  =  1.99  8  in  Ref.  20.   The  diameter  tolerance  allowed 
s 

at  the  throat  was  also  incurred  at  the  exit  plane  of  the 
nozzle.   The  area  ratios  of  the  nozzles  constructed,  there- 
fore, varied  between  1.6  82  3  and  1.69  79.   This  causes  a  Mach 

number  variation  between  M  =  1.996  and  M  =  2.001  for  one 

s  s 

dimensional,  inviscid  isentropic  flow  calculations,. 
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The  higher  Mach  number  nozzle  was  constructed  to  a 

nominal  area  ratio  of  2.6569  corresponding  as  above  to 

M  =  2.508.   Tolerances  in  manufacture  caused  a  possible 

spread  in  area  ratio  between  2.6366  and  2.6777  or,  in 

Mach  number,  between  M   =  2.500  and  2.516. 

s 

The  nozzles,  however,  are  actually  three  dimensional. 
Additionally,  the  viscosity  of  the  secondary  gas  will 
cause  boundary  layer  effects  to  be  present,  especially 
downstream  of  the  nozzle  throat  where  a  turbulent  boundary 
layer  will  exist.   The  nozzle  in  effect  will  sense  a 
smaller  area  ratio  than  that  for  which  the  nozzles  were 
designed.   The  overall  effect  is  to  decrease  Mach  number 
observed  from  that  determined  by  one  dimensional,  isen- 
tropic,  inviscid  calculations. 

Calibration  was  accomplished  by  suspending  a  thin  2.5° 
half-angle  cone  (with  its  apex  pointed  upstream)  in  the 
flow  exiting  the  nozzles.   The  apex  of  the  cone  was 
centered   in  the  exit  plane  of  the  nozzle.   Schlieren 
photographs  were  then  taken  of  the  flow  over  a  range  of 
nozzle  flow  stagnation-pressure-to-atmospheric-pressure 
ratios. 

The  conical  shock  angles,  6,  were  measured  from  the 
photographs  on  an  optical  comparitor.   The  angles  thus 
determined  were  6=  35.5°  ±  0.5°  for  the  low  Mach  number 
supersonic  nozzle  and  0  =  28.0°  ±  0.5°  for  the  higher  Mach 
number  supersonic  nozzle.   These  correspond,  respectively, 
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to  M  =  1.73  ±  0.05  and  M  =  2.13  ±  0.05  from  the  conical 
s  s 

shock  curves  of  Ref.  20. 

Taking  the  actual  uncertainty  in  M   caused  by  variation 
of  the  physical  dimensions  in  the  nozzles  to  be  of  the  same 
magnitude  as  that  determined  by  one-dimensional  flow 
calculations,  the  nozzles  constructed  operated  at  M  =1.73 
±  .002  and  2.13  ±  .008.   In  comparison  with  the  degree  of 
calibration  uncertainty  present,  the  dimensional  variation 
allowed  in  nozzle  construction  is,  therefore,  insignificant. 
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